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At Run II of the Tevatron it will be possible to measure the W boson mass with a relative precision of about 2 x 10 -4 , 
which will eventually represent the best measured observable beyond the input parameters of the SM. Proper interpretation 
of such an ultrahigh precision measurement, either within the SM or beyond, requires the meticulous implementation and 
control of higher order radiative corrections. The FORTRAN package GAPP, described here, is specifically designed to meet 
this need and to ensure the highest possible degrees of accuracy, reliability, adaptability, and efficiency. 



1. PRECISION TESTS 

Precision analysis of electroweak interactions follows 
three major objectives: high precision tests of the SM; 
the determination of its fundamental parameters; and 
studies of indications and constraints of possible new 
physics beyond the SM, such as supersymmetry or 
new gauge bosons. Currently, the experimental in- 
formation comes from the very high precision Z bo- 
son measurements at LEP 1 and the SLC, direct mass 
measurements and constraints from the Tevatron and 
LEP 2, and low energy precision experiments, such as 
in atomic parity violation, v scattering, and rare de- 
cays. These measurements are compared with the pre- 
dictions of the SM and its extensions. The level of 
precision is generally very high. Besides the need for 
high-order loop calculations, it is important to utilize 
efficient renormalization schemes and scales to ensure 
sufficient convergence of the perturbative expansions. 

The tasks involved called for the creation of a 
special purpose FORTRAN package, GAPP, short 
for the Global Analysis of Particle Properties [ [lj]. 
It is mainly devoted to the calculation of pseudo- 
observables, i.e., observables appropriately idealized 
from the experimental reality. The reduction of raw 
data to pseudo-observables is performed by the ex- 
perimenters with available packages (e.g., ZFITTER 
for Z pole physics). For cross section and asymme- 
try measurements at LEP 2 (not implemented in the 
current version, GAPP_99.7), however, this reduction 
is not optimal and convoluted expressions should be 
used instead. GAPP attempts to gather all available 
theoretical and experimental information; it allows the 
addition of extra parameters describing new physics; 
it treats all relevant SM inputs as global fit parame- 
ters; and it can easily be updated with new calcula- 
tions, data, observables, or fit parameters. For clarity 
and speed it avoids numerical integrations throughout. 



It is based on the modified minimal subtraction (MS) 
scheme which demonstrably avoids large expansion co- 
efficients. 

GAPP is endowed with the option to constrain non- 
standard contributions to the oblique parameters de- 
fined to affect only the gauge boson self-energies [ |) 
(e.g. S, T, and U); specific anomalous Z couplings; the 
number of active neutrinos (with standard couplings 
to the Z boson); and the masses, mixings, and cou- 
pling strengths of extra Z bosons appearing in models 
of new physics. With view on the importance of su- 
persymmetric extensions of the SM on one hand, and 
upcoming experiments on the other, I also included 
the b — > 57 transition amplitude, and intend to add 
the muon anomalous magnetic moment. In the latter 
case, there are theoretical uncertainties from hadronic 
contributions which are partially correlated with the 
renormalization group (RG) evolutions of the QED 
coupling and the weak mixing angle. These correla- 
tions will be partially taken into account by including 
heavy quark effects in analytical form; see Ref. [ | for 
a first step in this direction. By comparing this scheme 
with more conventional ones, it will also be possible to 
isolate a QCD sum rule and to rigorously determine 
the charm and bottom quark MS masses, m c and rhb, 
with high precision. 



2. GAPP 

2.1. Basic structure 

In the default running mode of the current version, 
GAPP_99.7, a fit is performed to 41 observables, out of 
which 26 are from Z pole measurements at LEP 1 and 
the SLC. The Fermi constant, Gf (from the muon life- 
time), the electromagnetic fine structure constant, a 
(from the quantum Hall effect), and the light fermion 
masses are treated as fixed inputs. The exception is 
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m c which strongly affects the RG runningF] of a{fi) for 
/i > rh c . I therefore treat fh c as a fit parameter and 
include an external constraint with an enhanced er- 
ror to absorb hadronic threshold uncertainties of other 
quark flavors, as well as theoretical uncertainties from 
the application of perturbative QCD at relatively low 
energies. Other fit parameters are the Z boson mass, 
Mz, the Higgs boson mass, Mjj, the top quark mass, 
rat, and the strong coupling constant, a s , so that there 
are 37 effective degrees of freedom. Given current pre- 
cisions, Mz may alternatively be treated as an addi- 
tional fixed input. 

The file f it . f basically contains a simple call to the 
minimization program MINUIT [ f§ (from the CERN 
program library) which is currently used in data driven 
mode (see smfit.dat). It in turn calls the core sub- 
routine fen and the x 2 -function chi2, both contained 
in chi2 . f . Subroutine fen defines constants and flags; 
initializes parts of the one-loop package FF [ || |(| ; and 
makes the final call to subroutine values in main.f 
which drives the output (written to file smf it . out). In 
chi2 the user actively changes and updates the data for 
the central values, errors, and correlation coefficients 
of the observables, and includes or excludes individual 
contributions to x 2 (right after the initialization, chi2 
= O.dO). To each observable (as defined at the begin- 
ning of chi2) corresponds an entry in each of the fields 
value, error, smval, and pull, containing the cen- 
tral observed value, the total (experimental and the- 
oretical) error, the calculated fit value, and the stan- 
dard deviation, respectively. The function chi2 also 
contains calls to various other subroutines where the 
actual observable calculations take place. These are 
detailed in the following subsections. 

Another entry to GAPP is provided through mh.f 
which computes the probability distribution function 
of Mjj- The probability distribution function is the 
quantity of interest within Bayesian data analysis (as 
opposed to point estimates frequently used in the con- 
text of classical methods) , and defined as the product 
of a prior density and the likelihood, £ ~ exp(— % 2 /2). 
If one chooses to disregard any further information 
on Mjj (such as from triviality considerations or di- 
rect searches) one needs a non-informative prior. It 
is recommanded to choose a flat prior in a variable 
defined on the whole real axis, which in the case of 
Mjj is achieved by an equidistant scan over log Mh- 
An informative prior is obtained by activating one of 
the approximate Higgs exclusion curves from LEP 2 
near the end of chi2.f . These curves affect values of 
Mjj even larger than the corresponding quoted 95% 
CL lower limit and includes an extrapolation to the 



kinematic limit; notice that this corresponds to a con- 
servative treatment of the upper Mjj limit. 

Contour plots can be obtained using the routine 
mncontours from MINUIT. For the cases this fails, 
some simpler and slower but more robust contour pro- 
grams are also included in GAPP, but these have to be 
adapted by the user to the case at hand. 

2.2. a, sm 2 9 w , M w 

At the core of present day electroweak analyses is 
the interdependence between Gf, Mz, the W boson 
mass, Mw, and the weak mixing angle, sin 2 9w- In 
the MS scheme it can be written as [ @, § , 
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s 2 is the MS mixing angle, c 2 = 1 — s 2 , and where, 



Af w = - A, 



and, 



n ww (M w ) -IW(0) 



V + B, (3) 



Afz = Afw + (1 — Afw 



nzz(Mf) 



Ml 



-, (4) 



collect the radiative cor- 

rections computed in sin2th.f. The II indicate MS 
subtracted self-energies, and V + B denote the ver- 
tex and box contributions to /i decay. Although these 
relations involve the MS gauge couplings they employ 
on-shell gauge boson masses, absorbing a large class of 
radiative corrections [ |i| . 

Afw and Afz are both dominated by the contribu- 
tion A~(Mz) which is familiar from the RG running 
of the electromagnetic coupling, 



a 



(5) 



1 Quantities defined in the MS scheme are denoted by a caret. 



and computed in alfahat.f up to four-loop O{ao i s ). 
Contributions from c and b quarks are calculated using 
an unsubtracted dispersion relation [ [| . If \x is equal to 
the mass of a quark, three-loop matching is performed 
and the definition of a changes accordingly. Pure QED 
effects are included up to next-to-leading order (NLO) 
while higher orders are negligible. Precise results can 
be obtained for \i < 2m n and fi > m c . 

Besides full one-loop electroweak corrections, Afw 
and Afz include enhanced two-loop contributions 
of 0(a 2 mf) [ [l(| (implemented using the ana- 
lytic expressions of Ref. [ |l)) and 0(a 2 m 2 ) \ |l2| 
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(available as expansions in small and large Mg); 
mixed electroweak/QCD corrections of 0(aa s ) [ |l3] ] 
and 0(aa 2 m 2 ) [ Q; the analogous mixed elec- 
troweak/QED corrections of 0(a 2 ); and fermion mass 
corrections also including the leading gluonic and pho- 
tonic corrections. 

2.3. Z decay widths and asymmetries 

The partial width for Z — ► // decays is given by, 

r // = %S^\Pf\ t 1 - ^\Qf\M^ 2 + 8Q^ 4 |%| 2 ] 



24s 2 c 



1 + 6, 



QED 



rNS 
°QCD 



.(6) 



P is the color factor, Q/ is the fermion charge, and 
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/5y and kf are form factors which differ from unity 
through one-loop electroweak corrections 
computed in rho.f and kappa. f, respectively. 



and are 
For 



f =/= b there are no corrections of 0(a mi) and con- 
tributions of 0{o?m\ ) to iif [ [l6| and pf [ [jj are very 
small and presently neglected. On the other hand, ver- 
tex corrections of 0(aa s ) [^) are important and shift 
the extracted a s by ~ 0.0007. 

The Z —> bb vertex receives extra corrections due 
to heavy top quark loops. They are large and have 
been implemented in bvertex.f based on Ref. [ |is|]. 
0{a 2 m\) corrections [ [n) are included, as well, 
while those of 0(a 2 m 2 ) are presently unknown. The 
leading QCD effects of 0(aa s m 2 ) and all sublead- 
ing 0(aa s ) corrections [ E(J are incorporated into pb 
and Kb, but not the 0(aoqm 2 ) contribution which is 
presently available only for nonsinglet diagrams [ . 

In Eq. (@), (5qed are the 0(a) and 0{a 2 ) QED cor- 
rections. <5qc D are the universal QCD corrections up 
to 0(a^) which include quark mass dependent contri- 
butions due to double-bubble type diagrams [ |2^, jiFJ] . 
()qcd are the singlet contributions to the axial-vector 
and vector partial widths which start, respectively, at 
0(a 2 ) and 0(al), and induce relatively large family 
universal but flavor non-universal mt effects [ Ej3l. 
The corrections appearing in the second line of Eq. (^|) 
are evaluated in leplOO.f. 

The dominant massless contribution to <5qcd can be 
obtained by analytical continuation of the Adler D- 
function, which (in the MS scheme) has a very well 
behaved perturbative expansion ~ 1 + J2i=o in 
a s = a s (Mz)/ir (see the Appendix for details). The 
process of analytical continuation from the Euclidean 
to the physical region induces further terms which are 
proportional to /3-function coefficients, enhanced by 
powers of ir 2 , and start at 0(al). Fortunately, these 
terms [ ^5) involve only known coefficients up to O(a^), 
and the only unknown coefficient in 0{ot & s ) is propor- 



tional to the four- loop Adler function coefficient, aV 
In the massless approximation, 



<rNS 
°QCD 



a s + 1.4092^ - (0.681 + 12.086)^ 
).19)a 4 s + {d A + 79.7)a 5 s + (d 5 - 121d. 



(7) 

3316)c4 
10 are clearly negligible. 



(g?3 — oa.iy )a s -|- ^114 -r iv.i )u s -|- [tb^ — ±4x113 

and terms of order a 7 s ~ 10" 
Notice, that the 0(oF s ) term effectively reduces the sen- 
sitivity to c?3 by about 18%. Eq. (0) amounts to a reor- 
ganization of the perturbative series in terms of the di 
times some function of a s ; a similar idea is routinely 
applied to the perturbative QCD contribution to r de- 
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of 0(m 2 f ) 



cays 

Final state fermion mass effects 
(and 0{mt) for b quarks) are best evaluated by ex- 
panding in m 2 (Mz) thus avoiding large logarithms 
in the quark masses. The singlet contribution of 
0{a 2 m 2 ) is also included. 

The dominant theoretical uncertainty in the Z line- 
shape determination of a s originates from the mass- 
less quark contribution, and amounts to about ±0.0004 
as estimated in the Appendix. There are several fur- 
ther uncertainties, all of O(10~ 4 ): from the 0(a 4 ) 
heavy top quark contribution to the axial-vector part 
of £q CD ; from the missing 0(aa 2 m 2 ) and 0(a 2 m 2 ) 
contributions to the Z6&-vertex; from further non- 
enhanced but cohering ©(aa^-vertex corrections; and 
from possible contributions of non-perturbative origin. 
The total theory uncertainty is therefore, 



Aa s (M z ) = ±0.0005, 



(8) 



which can be neglected compared to the current ex- 
perimental error. If rhb is kept fixed in a fit, then its 
parametric error would add an uncertainty of ±0.0002, 
but this would not change the total uncertainty (||). 

Polarization asymmetries are (in some cases up to 
a trivial factor 3/4 or a sign) given by the asymmetry 
parameters, 



l-4|g / |Rc(%)g 2 



1 -4\Q f \Re(k f )s 2 + . 



and the forward-backward asymmetries by, 
3 



AfbU) = -A e A 



f- 



(9) 



(10) 



The hadronic charge asymmetry, Qfb, is the linear 
combination, 



Q 



FB 



E 

q— e£,s,6 



R q A F B(q) 



R q A F B(q), 



(11) 



and the hadronic peak cross section, Chadj is stored in 
sigmah, and defined by, 

127rr e + e -r h ad 



°liad 



Mfr| 



(12) 
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Widths and asymmetries are stored in the fields 
gamma (f), alr(f), and afb(f). The fermion index, 
f, and the partial width ratios, R(f), are defined in 
Table |. 



Table 1 

Some of the variables used in leplOO . f . rj nv and I\ a d 
are the invisible and hadronic decays widths, respec- 
tively. 






V 


gamma (0) = r inv 


alr(0) 


1 


1 


e 


R(D = r had /r e+e - 


alr(l) = 


A e 


2 


M 


r(2) = ri ia d/r^+ M - 


alr(2) = 


A, 


3 


r 


R(3) = r ha d/r r+r - 


alr(3) = 


A T 


4 


u 


r(4) = rwr had 






5 


c 


r(5) = r ce /r ha d 


alr(5) = 


A c 


6 


t 


R(6) = 






7 


d 


r(7) = r dJ /r had 






8 


s 


r(8) = r sS /ri ia d 


alr(8) = 


A s 


9 


b 


RO) = r 65 /r h ad 


alr(9) = 


A b 


10 


had 


gamma(lO) = r had 


afb(10) = 


= Qfb 


11 


all 


gamma(ll) = Tz 







2.4. Fermion masses 

I use MS masses as far as QCD is concerned, but re- 
tain on-shell masses for QED since renormalon effects 
are unimportant in this case. This results in a hybrid 
definition for quarks. Accordingly, the RG running of 
the masses to scales ^ rh q uses pure QCD anoma- 
lous dimensions. The running masses correspond to 
the functions msrun(/i) , mcrunC/i), etc. which are 
calculated in masses. f to three-loop order. Anoma- 
lous dimensions are also available at four-loop order [ 
p8| , but can safely be neglected. Also needed is the 
RG evolution of a s which is implemented to four-loop 
precision [ |2^] in alfas.f. 

I avoid pole masses for the five light quarks through- 
out. Due to renormalon effects, these can be deter- 
mined only up to O(Aqcd) and would therefore in- 
duce an irreducible uncertainty of about 0.5 GeV. In 
fact, perturbative expansions involving the pole mass 
show unsatisfactory convergence. In contrast, the MS 
mass is a short distance mass which can, in principle, 
be determined to arbitrary precision, and perturbative 
expansions are well behaved with coefficients of order 
unity (times group theoretical factors which grow only 
geometrically). Note, however, that the coefficients of 
expansions involving large powers of the mass, iff 1 , are 
rather expected to be of 0(n). This applies, e.g., to 
decays of heavy quarks (n — 5) and to higher orders in 
light quark mass expansions. 

The top quark pole mass enters the analysis 



when the results on m t from on-shell produced top 
quarks at the Tevatron are included. In subroutine 
polemasses (nf ,mpole) m q (m q ) is converted to the 
quark pole mass, mpole, using the two-loop pertur- 
bative relation from Ref. [ ||(|. The exact three-loop 
result [ [jl| has been approximated (for mt) by employ- 
ing the BLM 



32 



scale for the conversion. Since the 
pole mass is involved it is not surprising that the coef- 
ficients are growing rapidly. The third order contribu- 
tion is 31%, 75%, and 145% of the second order for m t 
(nf = 6), nib (nf = 5), and m c (nf = 4), respectively. I 
take the three-loop contribution to the top quark pole 
mass of about 0.5 GeV as the theoretical uncertainty, 
but this is currently negligible relative to the experi- 
mental error. At a high energy lepton collider it will 
be possible to extract the MS top quark mass directly 
and to abandon quark pole masses altogether. 

2.5. v scattering 

The ratios of neutral-to-charged current cross sec- 
tions, 



T NC 
vN 
T CC ' 
'vN 



r NC 
vN 

T cc ■■ 



(13) 



have been measured precisely in deep inelastic v 
(9) hadron scattering (DIS) at CERN (CDHS and 
CHARM) and Fermilab (CCFR). The most precise re- 
sult was obtained by the NuTeV Collaboration at Fer- 
milab who determined the Paschos-Wolfenstein ratio, 



R~ = 



NC 
■N 



T cc 

vN 



T CC 

uN 



R v - ri? p , 



(14) 



with r = 



Io v n- 



' ' Results on R v are frequently 
quoted in terms of the on-shell weak mixing angle (or 
Mw) as this incidentally gives a fair description of the 
dependences on mt and Mh- One can write approxi- 
mately, 



where, 



gi + g A R- 

r 



= 9l-9w (15) 



2 1 



sin 2 6w + — sin 4 6w , g 2 R = -sa\ i 9w- (16) 
9 9 



However, the study of new physics requires the im- 
plementation of the actual linear combinations of ef- 
fective four-Fermi operator coefficients, £l,b.( u ) an d 
e L,R(d), which have been measured. With the ap- 
propriate value for the average momentum trans- 
fer, q2, as input, these are computed in the sub- 
routines nuh(q2,epsu_L,epsd_L,epsu_R,epsd_R) (ac- 
cording to Ref. [ |33|), nuhnutev, nuhccfr, and 
nuhedhs, all contained in file dis.f. Note, that the 
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CHARM results have been adjusted to CDHS condi- 
tions [ [54|. While the experimental correlations be- 
tween the various DIS experiments are believed to 
be negligible, large correlations are introduced by the 
physics model through charm mass threshold effects, 
quark sea effects, radiative corrections, etc. I con- 
structed the matrix of correlation coefficients using the 
analysis in Ref. [ pil , 





Ru 


Ru 


Ru 


Re 


R r , 


Rd 


/ 1.00 


0.10 


0.10 


0.10 


0.00 


0.00 


0.00 \ 


0.10 


1.00 


0.40 


0.40 


0.10 


0.10 


0.10 


0.10 


0.40 


1.00 


0.40 


0.10 


0.10 


0.10 


0.10 


0.40 


0.40 


1.00 


0.10 


0.10 


0.10 


0.00 


0.10 


0.10 


0.10 


1.00 


0.15 


0.15 


0.00 


0.10 


0.10 


0.10 


0.15 


1.00 


0.15 


V 0.00 


0.10 


0.10 


0.10 


0.15 


0.15 


1.00 j 



.(17) 



The effective vector and axial-vector couplings, gy e 
and <7^ e , from elastic ve scattering are calculated in 
subroutine nue(q2,gvnue,ganue) in file nue.f. The 
momentum transfer, q2, is currently set to zero [ 
j36| . Needed is the low energy p parameter, rhonc, 
which describes radiative corrections to the neutral- 
to-charged current interaction strengths. Together 
with sin2t0 (described below) it is computed in file 
lowenergy . f . 

2.6. Low energy observables 

The weak atomic charge, Qw, from atomic parity 
violation and fixed target ep scattering is computed 
in subroutine apv(Qw,Z,AA,Clu,Cld,C2u,C2d) where 
Z and AA are, respectively, the atomic number and 
weight. Also returned are the coefficients from lepton- 
quark effective four-Fermi interactions which are cal- 
culated according to [ [57) . 

These observables are sensitive to the low energy 
mixing angle, sin2t0, which defines the electroweak 
counterpart to the fine structure constant and is sim- 
ilar to the one introduced in Ref. [ M. There is sig- 
nificant correlation between the hadronic uncertainties 
from the RG evolutions of a and the weak mixing an- 
gle. Presently, this correlation is ignored, but with the 
recent progress in atomic parity violation experiments 
it should be accounted for in the future. 

An additional source of hadronic uncertainty is in- 
troduced by 7-Z-box diagrams which are unsuppressed 
at low energies. At present, this uncertainty can be 
neglected relative to the experimental precision. 

Besides apv, file pne . f contains in addition the sub- 
routine moller for the anticipated polarized fixed tar- 
get M0ller scattering experiment at SLAC. Radiative 
corrections are included following Ref. [ B8| . 



2.7. b -> S7 

Subroutine bsgamma returns the decay ratio, 

B(b -> 8j) 
B{b cev) ' 

It is given by [ ||, E9], 



R 



6 a 



Vcl 



S \D\ 2 +A/S + Snp + Sew 

W) (i + ^p)(i + ^) ■ 



(18) 



(19) 



where \V^ S Vtb/V c b\ 2 = 0.950 is a combination of 
Cabbibo-Kobayashi-Maskawa matrix elements and S 
is the Sudakov factor [ ^l). Snp and Sew ar e non- 
perturbative and NLO electroweak [ ^] corrections, 
both for the b — > sj and the semileptonic (b — ► cev) 
decay rates. 



D 



a a (m b ) tril 
°7 + 4?r i°7 



n 



(20) 



is called the reduced amplitude for the process b — > sj, 
and is given in terms of the Wilson coefficient C 7 at 
NLO. C 7 and the other d appearing below are ef- 
fective Wilson coefficients with NLO RG evolution [ 
fill from the weak scale to fx = rhb understood. The 
NLO matching conditions at the weak scale have been 
calculated in Ref. [ Q . D includes the virtual gluon 
corrections, 

V = r 2 C 2 ° + r 7 C 7 °+r 8 C 8 , (21) 

so that it squares to a positive definite branching frac- 
tion. On the other hand, the amplitude for gluon 
Bremsstrahlung (b — > sjg), 

A = ^^[C 2 (C 8 °/28(1) + C 7 °/ 27 (l) + C 2 °/ 22 (l))+ 

C 8 °(C 8 / 88 (<5) + C 7 °/ 7 8(^)) + (C 7 °) 2 /77(<5)], (22) 

is added linearly to the cross section. The Wilson co- 
efficient C 2 is defined as in Ref. [ E5J . It enters only 
at NLO, is significantly larger than C^, and dominates 
the NLO contributions. The parameter < S < 1 
in the coefficient functions characterizes the min- 
imum photon energy and has been set to 5 = 0.9 [ 
EjJ , except for the first line in Eq. (22) where S = 1.0 
corresponding to the full cross section. The / 2 j are 
complicated integrals which can be solved in terms of 
polylogarithms up to 5th order. In the code I use an 
expansion in z = m^/m^ and S — 1.0. Once experi- 
ments become more precise the correction to 5 = 0.9 
should be included. 

f(z) is the phase space factor for the semileptonic 
decay rate including NLO corrections [ fl6|. I defined 
the MS mass ratio in z = [m c (mb)/?7i{,(mf,)] 2 at the 
common scale, fi = rhb, which I also assumed for the 
factor ml multiplying the decay widths. Since I do not 
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rccxpand the denominator this effects the phase space 
function at higher orders. Using the 0(a 2 )-estiniatcP] 
from Ref. [ [47|, I obtain for the semileptonic decay 
width, 



(23) 



where fo(z) is the leading order phase space factor. It 
is amusing that the coefficients in Eq. ( p3| ) are comfort- 
ably (and perhaps somewhat fortuitously) small, with 
the O(a^)-coefficient even smaller than the one in Ref. [ 
|f7[ where a low scale running mass had been advo- 
cated. Moreover, using the prefactor 7tt| in the numer- 
ator of R reduces the size of r-j in Eq. and there- 
fore the coefficient n(6) = /77(d) + rj/2 which multi- 
plies the term a s (C°' cff ) 2 . I obtain -2.1 < k(S) < 1.4, 
while with the pole mass prefactor Mjj one would have 
-8.7 < k(S) < -5.3. 
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A. Uncertainties from perturbative QCD 

Writing the perturbative expansion of some quantity 
in its general form for an arbitrary gauge group, it can 
easily be decomposed into separately gauge invariant 
parts. Table || shows for some (related) examples that 
after removing the group theoretical prefactors, all co- 
efficients, i/i, are strictly of order unity, and that their 
mean, y, is very close to zero. In particular, there is no 
sign of factorial growth of coefficients. These observa- 
tions offer a valuable tool to estimate the uncertainties 
associated with the truncation of the loop expansion, 
so I would like to make them more precise. 

Assume (for simplicity) that the yi are random draws 
from some normal distribution with unknown mean, 
fx, and variance, a 2 . One can show that the marginal 
distribution of \x follows a Student-t distribution with 
n — 1 degrees of freedom, t n —i, centered about y, and 
with standard deviation, 



T«(yi-i/) 2 

n(n — 3) 



(24) 



As can be seen from the Table, \i is consistent with zero 
in all cases, justifying the nullification of the unknown 
coefficients from higher loops. I next assert that the 



2 1 computed the 0(a 2 ) coefficient for comparison only, and did 
not include it in the code. 



distribution of a, conditional on /j = 0, follows a scaled 
inverse-^ 2 distribution with n degrees of freedom, from 
which I obtain the estimate, 



a = erg i Act 



i± 



1 



2 (71 - 4) 



(25) 



Inspection of the Table shows indeed that a, as the 
typical size of a coefficient, is estimated to be < 0(1). 

I now focus on the par tial hadronic Z decay width. 
As discussed in Section |2.3| , the O(o%) term, cfe, is 
much smaller than the 7r 2 term arising from analytical 
continuation. This is specifically true for the relevant 
case of n/ = 5 active flavors, where large cancellations 
occur between gluonic and fermionic loops. Notice, 
that the D- function, in contrast to -Rhacb has opposite 
signs in the leading terms proportional to C A Cp and 
CaCfTfUf. Indeed, the Adler D-function and the 
/3-function have similar structures regarding the signs 
and sizes of the various terms (see Table ||), and we 
do expect large cancellations in the /3-function. The 
reason is that it has to vanish identically in the case of 
N = 4 supersymmetry. Ignoring scalar contributions 
this case can be mimicked by setting TpUf = 2Ca 
(there are 2 Dirac fermions in the N — 4 gauge multi- 
plet) or n/ = 12 for QCD, which is of the right order. 
In fact, all known QCD /3-function coefficients become 
very small for some value of 71/ between 6 and 16. We 
therefore have a reason to expect that similar cancel- 
lations will reoccur in the di at higher orders. As a la 
error estimate for d 3 , 1 suggest to use the largest known 
coefficient (3 x 0.71) times the largest group theoretical 
prefactor in the next order (C a Cf) which results in 



0± 77. 



(26) 



With Eq. (0) and a s (M z ) = 0.120 one can absorb all 
higher order effects into the ©(a^-coefficient of i?i ia d, 
r§ ff = —81 ± 63. This shifts the extracted a s from 
the Z lineshape by +0.0005 and introduces the small 
uncertainty of ±0.0004. 

The argument given above does certainly not apply 
to the quenched case, n/ = 0, and indeed diirif = 0) is 
about —73% of the 7r 2 term, i.e., large and positive. In 
the case of 71/ = 3, which is of interest for the precision 
determination of a s from t decays, di is about —38% 
of the 7r 2 term. If one assumes that the same is true 
of d 3 , one would obtain d^(nf = 3) = 60. Estimates 
based on the principles of minimal sensitivity, PMS, 
or fastest apparent convergence, FAC, yield d^ijif = 
3) = 27.5 [ |2q| so there might be some indications for a 
positive ds(nf = 3). In any case, all these estimates lie 
within the uncertainty in Eq. ( |26| ) and we will have to 
await the proper calculation of the 0(aj)-coefficient to 
test these hypotheses. Note, that the current r decay 
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analysis by the ALEPH Collaboration uses — 50 ± 
50 [ |5(| which is more optimistic. 

The analogous error estimate for the five-loop (3- 
function coefficient yields, 



th 



0± 579. 



(27) 



To get an estimate for the uncertainty in the RG run- 
ning of & s , I translate Eq. (|27|) into 



P3 — P3 ± PA, 



(28) 



where fia is taken to be the lowest scale in- 
volved. This overestimates the uncertainty from /?4, 
thereby compensating for other neglected terms of 
0(a™ +4 In" £i 2 /^o)- For the RG evolution from fx — m T 
to fi = Mz this yields an uncertainty of Aa s (Mz) = 
±0.0005. Conversely, for fixed a s {M z ) = 0.120, I ob- 
tain a s (rh b ) = 0.2313 ± 0.0006, a s (m T ) = 0.3355 ± 

0. 0045. and d s (m c ) = 0.403±0.011, where I have used 
rhf, = 4.24 GeV and m c = 1.31 GeV. For comparison, 
the ALEPH Collaboration quotes an evolution error of 
Aa s (Mz) — ±0.0010 which is twice as large. I empha- 
size that it is important to adhere to consistent stan- 
dards when errors are estimated. This is especially true 
in the context of a global analysis where the precisions 
of the observables enter as their relative weights. 
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Table 2 

Coefficients (MS ) appearing in the /3-function of a 
simple group [ ^9|; in non-Abelian corrections to the 
QED /3-function (denoted D) [ |, ||; in the Adler 
D- function [ |^| (rescaled by an overall factor 1/3); 
and in i?h a( j (analytical continuation of D). The first 
four segments correspond, respectively, to the first four 
loop orders of non-singlet type. The fifth segment is 
the singlet (double triangle) contribution in 0{o? s ). In 
D, D, and i?had, an overall factor aMz and the sums 
involving charges or Z couplings have been dropped. 
The completely symmetrical tensors of rank four, dA 
and d F , as well as Tp when appearing in parenthesis, 
apply to the /3-function only. Each Tp is understood to 
be multiplied by the number of flavors, n/, except for 
the singlet term involving the symmetrical structure 
constants, d. 



group factor 





D 


D 


-Rhad 


Ca 


0.92 








\ 1 F) 


yj.oo 


— n 

u.oo 


u.oo 


u.oo 


rf 1 


0.71 








CaTf 


-0.42 








C F {Tp) 


-0.25 


-0.25 


0.25 


0.25 


c \ 


0.83 








C\T F 


-0.82 








C'aCf(T f ) 
C'aT f 


-0.36 


-0.23 


0.18 


0.18 


0.09 










0.03 


0.03 


-0.03 


-0.03 


C F T F {Tp) 


0.08 


0.08 


-0.06 


-0.06 


C% 


1.19 








C\T F 


-1.67 








C A C F (T F ) 
r 2 r 2 

L/ A 1 F 

CaC f (T f ) 


-0.23 


-0.28 


0.32 


-0.38 


0.50 
-0.42 


-0.07 


0.51 


0.51 


C'aC f T f (T f ) 
C A T F 


0.51 
0.01 


0.42 


-0.71 


-0.21 


C 3 f (Tf) 

C 2 F T F {T F ) 

C F T F (T F ) 


0.18 


0.18 


-0.18 


-0.18 


-0.17 


-0.17 


0.02 


0.02 


0.02 


0.02 


0.09 


-0.01 


d 2 A /N A 
d A d F /N A 


1.07 
-2.38 








d F /N A (T F d 2 /4) 


0.50 


0.50 


-0.50 


-0.50 


y 


-0.02 


-0.01 


0.02 


-0.01 


A/x 


0.17 


0.09 


0.11 


0.09 


cro 


0.83 


0.28 


0.37 


0.31 


Act 


0.13 


0.07 


0.09 


0.08 
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